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The characterization of the transverse structure of the QCD string is discussed. We formu- 
late a conjecture as to how the stress-energy tensor of the underlying gauge theory couples 
to the string degrees of freedom. A consequence of the conjecture is that the energy density 
and the longitudinal-stress operators measure the distribution of the transverse position of 
the string, to leading order in the string fluctuations, whereas the transverse-stress operator 
does not. We interpret recent numerical measurements of the transverse size of the confining 
string and show that the difference of the energy and longitudinal-stress operators is the 
appropriate probe to use when comparing with the next-to-leading order string prediction. 

Secondly we derive the constraints imposed by open-closed string duality on the transverse 
structure of the string. We show that a total of three independent 'gravitational' form factors 
characterize the transverse profile of the closed string, and obtain the interpretation of recent 
effective string theory calculations: the square radius of a closed string of length /? defined 
from the slope of its gravitational form factor, is given by log in d space dimensions. 
This is to be compared with the well-known result that the width of the open-string at 
mid-point grows as log ^ . We also obtain predictions for transition form factors among 
closed-string states. 

PACS numbers: 11. 15. Ha, 11.25.Pm, 12. 38. Aw 

I. INTRODUCTION 

The area law of Wilson loops in lattice gauge theories [TJ has long been interpreted in terms 
of a string formation by the flux lines. In SU(3) gauge theory the law, (W) ~ e~ uA , si gnals 
the linear confinement of heavy quarks Q and Q: the static potential takes the form V(r) ~ err, 
where a is identified with the string tension. Once the quark-interdistance r is significantly larger 
than the confinement scale -y/cx, it was realized a long time ago that the corrections to the static 
potential, as well as the low-energy excitations of the QQ system, could be described by an effective 
two-dimensional theory [2]. This 'worldsheet' theory of the d — 1 [27] massless degrees of freedom 
h, namely the transverse fluctuations of the string, led to two important predictions: firstly, the 
linear potential receives a 1/r corrections, the Liischer term [3], and its excitations are spaced by 
- gaps. Secondly, the amplitude of the transverse string fluctuations grows logarithmically with 
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the length of the string [4], 

Wo-(^) = ^lo g -. (1) 

lira ro 

It is this second aspect of the low-energy string dynamics that is the focus of this paper. 

Recently, highly accurate numerical results have been obtained in the d = 2 SU(2) gauge theory 
for the expectation value of local operators in the presence of a static QQ pair [5]. Measured as a 
function of the distance \y\ from the QQ axis, is defines a distribution whose second moment was 
successfully compared to the effective-theory prediction ([!]). In view of these results and of the 
prospect of pushing the comparison to next-to-leading order, the first issue we wish to address is the 
precise connection between the profile probed by a local gauge-theory operator and the worldsheet 
expectation value of h 2 . 



II. COUPLING OF THE STRESS-ENERGY TENSOR TO THE CONFINING STRING 

With the physical picture of a fluctuating 'thin' string in mind, the stress and energy stored 
in the flux lines is entirely carried by the string. If only the transverse component of a string 
element's motion contributes to the string energy, then the Hamiltonian is H = f , dm , where 
dm is the rest mass of an element of the string. For the Nambu-Goto string, dm = ads, where ds 
is the length of the string element, but more sophisticated possibilities, such as a curvature term, 
should be kept in mind [6j[28 . Retaining only the simplest rest mass contribution, the expression 
for the string energy density in Minkowski space reads, in the static gauge, 

T 00 (t, yi ,y) = a - l + {dih) 2 d-1 _ 

Here h is the worldsheet field; it has d — 1 components. The worldsheet is parametrized by 
t and j/i, the coordinate that runs along the QQ axis, and y contains the d — 1 coordinates 
transverse to the string. The worldsheet indices are denoted generically by a, b, . . . We thus expect 
the energy-density operator of the underlying gauge theory to couple to the worldsheet operator 
appearing in this expression, understood as an expansion in d a h. At zeroth order, we have simply 
Too(t,yi,y) = a5 d ^ 1 (y — h(t,yi)), which means that the distribution measured by the energy 
density operator coincides with the distribution in h. In particular the second moments in y of the 
transverse distribution obtained from Too is expected to match the expression for (h 2 ) calculated 
in the worldsheet theory. Expanding Eq. ([2]), we obtain 

T 00 = a8(y-h)\l^((^ t h) 2 ^^ 1 h) 2 )-l[(^ 1 h) 2 } 2 ^[(^^ (3) 
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This expression suggests that at leading order in the fluctuations, J d d 1 y y 2 Too(t, yi,y) measures 
the worldsheet expectation value of the operator 

h^l + Udthf + l&h) 2 ). (4) 

Thus when comparing Monte-Carlo data for / d d ^y y 2 T 00 (t, yi , y) with the effective theory, the 
worldsheet expectation value of h 2 needs to computed to next-to-leading order, a tour-de-force 
achieved very recently [7], but also the leading order expectation value of h 2 ^((dth) 2 + (d\h) 2 ) 
needs to be calculated. It is probably simpler to work with the operator Too — f° r which we 



will see that the undesirable contribution of the quadratic fluctuations cancels out (Eq. 12). 

It is instructive to note that the energy density expression ^ derived from geometric consid- 
erations coincides with the form of the canonical energy density derived from the Luscher-Weisz 
worldsheet action with the standard Noether procedure. Indeed the NLO Lagrangian reads 

C ws = \d c h ■ d c h + c 2 (d a h ■ d a h){d b h ■ d b h) + c 3 {d a h ■ d b h){d a h ■ d b h) + ... (5) 

with a priori free coefficients c 2 and C3, and the stress-energy tensor 

rj = d a h ■ d b h + A C2 {d c h ■ d c h) {d a h ■ d b h) + 4 C3 (d a fc • d c h) {d b h ■ d c h) - g ab c ws , (6) 

with in particular 

= |((5 i /i) 2 +(9 1 /i) 2 )+(c 2 + C 3)((a 1 /i) 2 -3(^/i) 2 )+2 C2 (a^) 2 (a 1 /i) 2 +2 C3 (ai/i-5 1 /i) 2 . (7) 

Expressions tlh and fclti are consistent for c 2 = I and C3 = —4, which are the Nambu-Goto 



8 ^ ~ 4' 

values 



In d = 2 space dimensions, it was shown in [8] that these values for the 'low-energy constants' 
c 2 and C3 are the only ones compatible with open-closed string duality. It was subsequently shown 
that this requirement is also equivalent to requiring that closed string have a relativistic dispersion 
relation, in other words requiring Poincare invariance [9|. If one requires that the effective string 
theory also describes a situation where the worldsheet itself is a torus in a way that is consistent with 
the open- and closed-string spectral representations, then these values are the only ones possible 
in any dimension [ID]. In view of the geometric interpretation of the energy density operator, 
these results show that only a string that is 'immaterial', i.e. for which only transverse motion 
of an element of the string contributes to the string energy, yields a spectrum that is consistent 
with open-closed string duality. Were it not for this fact, the fraction in Eq. Q would have been 
replaced by \J i^q^ z, which in particular does not yield mixed term (dth • d\h) 2 . In other words, 
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numerical evidence that the open string spectrum requires c 2 and C3 to take up their respective 
Nambu-Goto values really confirms the 'immaterial' nature of the confining string. 

It is also of interest to write out the expressions for the longitudinal stress operator Tn explicitly 
(see Eq. @), 

T n (t, yi ,y) = a5(y - h(t, yi )) [ -1 + \{{d t h) 2 + {dthf) + (c 2 + c^^h) 2 - (d t h) 2 ) 

-2c 2 (d^) 2 (difr) 2 " 2c 3 (d t h ■ di/1) 2 ] . (8) 

This formula implies that the transverse string profiles obtained with Too and Tn differ at quadratic 
order in h. There is a specific reason why Tn is an interesting probe of the string profile. The 
transverse profile of the open string depends in general at what point y\ along the string it is 
measured. It is easy to see that if one uses Tn, then the total longitudinal stress inside a transverse 
spatial slice, f d d ~ 1 yTn(yi,y), does not depend on the position y\ of the slice along the string. 
This is simply because from the closed-string point of view, Tn plays the role of the energy density 
operator, and therefore its forward matrix elements are diagonal in an energy-eigenstate basis. 
Evaluated on the ground state, this integrated longitudinal stress yields the static force, 

J d d - 1 y(T 11 (y 1 ,y)) QQ = -^l (V yi ). (9) 

Because of this distinguishing property, the integrated longitudinal stress is conserved along the 
open string, and it is natural to ask how the transverse distribution of longitudinal stress changes 
as one moves along the string. 

It is however not true than any operator tracks the movement of the string at leading order. Take 
for instance the transverse operator T22. There is no corresponding operator on the worldsheet, 
since it is a two-dimensional field theory. One can show that 

J d m J d d - 1 y(T 22 (y 1 ,y)) QQ = 0, (10) 

when correlated with the pair of Polyakov loops. The physical reason why the three-point function 
of T22 vanishes is that the string does not, on average, exert any stress along the transverse 
directions. Because the sum rule of this operator does not yield a term proportional to the length 
of the string, this operator is not measuring, to leading order in the fluctuations h, the position of 
the string. Therefore one cannot define a transverse distribution of the string with a probabilistic 
interpretation based on this operator. Instead this operator is sensitive in leading order to the 
expectation value of higher-derivative worldsheet operators. 

We have followed the approach of Liischer and Weisz [8j and worked in the static gauge. The 
point of view adopted by Polchinski and Strominger |llj puts more emphasis on the conformal 
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symmetry of the worldsheet theory, which severely constrains the class of actions they consider. 
It is therefore worthwhile to investigate the fate of conformal symmetry in the static gauge as 
well. This issue is left for a future study. We simply note that the trace of the canonical energy- 
momentum tensor 

T™> a = 2c 2 (d a h ■ d a hf + 2c 3 {d a h ■ d b h)(d a h ■ d b h) + ■■■ = 2£< 4 > + . . . (11) 

no longer vanishes at the quartic order. However it is well-known that the canonical energy- 
momentum tensor is in general not traceless even when the field theory is conformally invariant. 
It can however be improved [12J in the sense that terms A a b that satisfy d a A a t> = and do not 
modify the conserved charges can be added in such a way that T ab is traceless when the theory is 
conformal. See [T3] for a discussion in two-dimensional field theory. It would be interesting to see 
whether the line of low-energy constants C3 = — 2c 2 |8j plays a special role in this respect. 
An observation of 'practical' importance is that the linear combination 

(Too - T n )(t, yi,y) = 2a5(y - h(t, y x )) (l + O(dh) 4 ) (12) 

tracks the position of the string up to possibly quartic corrections. This makes it the simplest 
operator to measure the mean square radius of the string fluctuations, which can be compared 
directly with the NLO formula obtain in [7j. 

The rest of this paper is structured as follows. We start by studying the structure of the 



confining string as seen by the energy-momentum tensor in section (III). We then work out the 
constraints on three-point correlation functions imposed by the open-closed string duality in section 
(IV). The leading-order string formula 0, generalized to contain the contributions of excited states 



in the three-point function, turns out to be consistent with the functional form in r imposed by the 
closed-string spectral representation, and we thereby identify the effective theory prediction for the 
form factors of the closed strings. In particular we find that the square radius of the ground state 
closed string, defined in the standard way from the slope of its form factors at the origin, grows 
logarithmically with the length of the string /3. In section |V| ) we give the explicit form of the 
energy-momentum tensor on the lattice in d + 1 dimensions. This allows us to interpret a recent 
high-accuracy calculation of the string width in numerical lattice gauge theory in terms of matrix 
elements of the energy-momentum tensor. In the rest of this paper, we work in Euclidean space, 
and our sign conventions are as follows. In Minkowski space, the thermal expectation values of the 
diagonal components are (Too) = e and (Tn) = p (respectively the energy density and pressure), 
while in Euclidean space (Too) = e and (Tn) = —p. 
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III. GRAVITATIONAL FORM FACTORS OF CLOSED STRINGS 

In this section, we analyze how the transverse size of closed strings can be characterized. In 
the pure SU(iV) gauge theory, the only conserved charges are energy and momentum. Therefore, 
it is natural to measure the width of the string in terms of the distribution of these charges. While 
for the open string, the width can be probed directly in x-space, it has to be defined initially in 
momentum space through a form factor for the closed string: the form factors with respect to the 
energy- momentum tensor T^ v are the Fourier transforms of the energy and distributions. This 
simple relation between form factors and charge distribution applies because of the non-relativistic 
kinematics of the closed string, by which we mean that their transverse size is parametrically larger 
than their inverse mass. By contrast, the electromagnetic form factors of the proton only correspond 
to the Fourier transform of charge and magnetization in the infinite- momentum frame |14j . 

Here we will restrict ourselves to studying the form factors of states that contribute to the 
Polyakov loop two-point function. These states are translationally invariant in the longitudinal 
direction, therefore we restrict the momentum transfer to the transverse directions. Furthermore, 
the closed string states are rotationally invariant, hence they have spin zero in (d— l)-dimensional 
space. 

In order to exhaustively list the relevant form factors, we decompose the full (d + l)-dimensional 
energy-momentum tensor into irreducible representations of (i-dimensional space. The closed 
strings are stretched around a cycle of length /3 in a spatial direction labeled z, while the other 
spatial directions are labeled by k, I = 1, . . . , (d — 1). Schematically, the decomposition takes the 
form 



Too Tot 
Tko T k i 



Tqz 



\ 



T, 



kz 



T 

± zz 



(13) 



\ T z o T z k 

In the following, we choose the normalization of states such that 



(tP, P'\i>, P) = {2^) d - l 5 d -\P' -P)-P- 2E P . (14) 

The operator T zz , which measures the stress in the z direction, is a scalar operator from the point 
of view of physics within an z = constant slice. Therefore, its matrix elements can be parametrized 
as 

<V/,P'|T zz (0)|V,P) = 2MM'h{^'^q 2 ). (15) 
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We use the standard notation q = P' — P, P = |(P + P'), and have accounted for the possibility 
that the mass of the final state M 1 differs from the mass of the initial state M. 

Secondly, we note that (Tq z , Tf. z ) is a conserved vector from the point of view of a z = constant 
slice, if one restricts oneself to matrix elements between states that are translationally invariant in 
the z direction: 

d T 0z + 8 k T kz + d z T zz = 0. (16) 
=o 

We are thus dealing with the vector form factor of a scalar object, hence (by analogy with the pion 
electromagnetic form factor), 

W,P'\T flz (0)\i;,P) = MP fl f(i;',i;;q 2 ), y,±z. (17) 

However, Tn Z is odd under the reflection z — > —z. For matrix elements with P z = P z = 0, this 
implies that / must vanish identically |30| . 

Finally, the components of T^ v not containing the index l z' form a tensor with respect to the 
SO(d) group. Taking again into account the fact that these components form a conserved tensor in 
the subspace of states invariant under translations along the y direction, one finds that the general 
form of the matrix elements of T^ u is 

{i//, P'lT^O)^, P) = -2jy\,/i(^, V; q 2 ) + 2{q^q v - q\ u ) f 2 (^ , ^ q 2 ), fi, v + z. (18) 

Thus the transverse structure of the ground state of the string is characterized by a total of three 
form factors {/i}| =1 - The matrix elements 

(V>,P|V(0)|V>,P> = -2P^P„ (19) 

dE 2 

(iP,P\T zz (0)\iP,P) = (20) 
determine the forward, diagonal matrix elements of /i and f%, 

/i(V,V;0) = i, /3(^;^;0) = 2 ^ 2 ^ - (21) 

The interpretation of these form factors is that measures the transverse distribution of longitu- 
dinal stress in the string, while f\ measures the transverse distribution of energy. The form factor 
$2 is somewhat less obvious to interpret. For two states with momenta equal and opposite aligned 
along the direction 1 (Breit frame), it describes the ability of T22 to induce a transition between 
these states per unit (momentum-transfer) 2 (this interpretation requires d > 3). Indeed, in this 



kinematic configuration, f% does not contribute to the matrix element ( 18 ) 
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y, x - y 









FIG. 1: The geometry of the three-point function. 



IV. TRANSVERSE STRUCTURE OF OPEN AND CLOSED STRINGS 

The goal of this section is to derive the spectral representation of a three-point function where 
a local operator is used to probe the structure of the confining string. We begin by recalling the 
spectral representation of the Polyakov loop two-point function. The geometry of the Polyakov 
correlator is illustrated in Fig. Q. 

The open-string representation of the Polyakov loop two-point function reads, setting r 2 = 
x\ + x 2 and with w n integer weights, 

(P o (x 1 ,x)P$(0,O)) =J2^ne- V ^ r)L . (22) 

n 

Upon introducing the matrix elements |31j 

b n = (vac|P (0,O)|n,p), (23) 
the closed-string representation of the same correlation function reads 

m^)PS M ) = -^\K?j J^^^y (24) 

„ yKt (MA^ e ~M n r. (26) 
In the last line we have used the asymptotic form of the modified Bessel function, K v (x) ~ e~ x — • 



2x ' 

the result is equivalent to using non-relativistic kinematics to begin with. This expression dictates 
the functional dependence on r of the Polyakov loop correlator. As usual in deriving relations 
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between open and closed strings, the correlation function cannot be simultaneously dominated by 
a single open-string state and a single closed-string state. Let j3 be the length of the closed strings. 
For (3 3> r, a single open-string state dominates, but 0(/3/r) closed-string states contribute in 



Eq. (|26|). 

Consider now the connected correlation function of a pair of Polyakov loops in the direction 
and a local operator O. Figure M illustrates the geometry of the correlator. Its spectral 



interpretation in terms of open-string states reads, for f3 S> a 2 . 



(P (x u x) O(y , yi ,y) P *(0, 0)) = £ e~ v ^ (0(y u y)) n (27) 



In terms of closed-string states it can also be written as 

1 d— 1 „/ r a <l— 1 , 



(P (x 1 ,x)O(y ,y 1 ,y)P^0,0)) = J ±^e^ / A^e^/lp', g, x l5 (28) 



where 



/(p',q,x 1 ,y 1 ) = I d^W^y d^y e^^^i.x) O(lft,,W,»)i^(0,0)) (29) 

is the correlation function in momentum space, which has a more natural interpretation from the 
closed-string point of view. Here p and q have d— 1 components. Due to the translation invariance 
of the Polyakov loops along the direction, / has no dependence on yo, which we therefore choose 
to be zero. 



With the normalization of states given by Eq. (14), we parametrize the matrix elements by 

(m,p'\0\n,p) = 2M m M n F m ' n (p,q), (30) 
p = p'-q, p = p' -\q=\{p + p'). (31) 

This parametrization is designed for dimension (d+1) operators, for which F m ' n is dimensionless. 
We can now write the spectral representation of /, 

/(f/,g,si,tt)=;jj> 2Em{pf)p 2M m M n ^ (p , q) -- mK , (p = p>-q). (32) 

Next we specialize to the case of a scalar operator with respect to the symmetry group SO(d) of a 
time-slice. Examples thereof are Too or T w . We will return to the case of an operator with a more 
general tensor structure in section ( |IVD[ ). Thus F is a function of q 2 alone, hence 

(P (x 1 ,x)O(y ,y 1 ,y)P^0,0)) = £ ^ f -^-^ e -*» 2M m M n F m < n (q 2 )l mn ( yi , x x - Vl , x, q), 

m,n " J \ ■> 

/(ifi-lpl e ~ip' '■x-Em{p')y'2-E n (p' '-q)yi 
(2^ 2E m (p')2E n (p / — q) " (33) 
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We choose without loss of generality x = 0. The quantity X is a massive one-loop integral, 



2^(2/1, 2/2, 0,q) 



oV 



d d-i p 



1 



1 



(34) 



2vr " J 2vr " 7 (2tt) < '- 1 w 2 + (p - g) 2 + M 2 w /2 + p 2 + M 2 , ' 
This integral can be treated by standard techniques of quantum field theory, see for instance [TS] 
p. 327. However we anticipate that non-relativistic kinematics is sufficient to study the long-distance 
behavior of the correlators (in the effective string theory, this will be guaranteed as long as <r/3yi 3> 
1), 

M m M n \ exp -(M m y 2 + M nVl + £ ^ 



i m n(yi,y2,o,q) ~ , , , , . 

\V2M n + yiM m 
Therefore, with y 2 = x\ — y±, 



2 M myi +M n y 2 - 



(2vr)5^ 1 ) 2M m -2M n 



b m bt 



<Po(xi,0)0(yo,yi,y)Po*(0,0)) ~ E ^ 



M m M„ 



2/ 2 M n + yi M„ 



p -(i/sM m +!/iM„) 



2(2vr) 



(35) 



(36) 



_ ?_p-iQ y F m ' n (a 2 }p~ 

( 27r )d-l e f 19 je 



WH/2 



2 M m! , 1 +M„i/2 



This expression dictates the leading-order functional dependence on x\ and y% of the three-point 
function that the effective string theory must respect. 



Expression (36) can be viewed as a distribution in y. The quantity we will confront with a 



prediction from the effective string theory is its second moment at x\ = r, 



w 2 (r,/3, yi ) = 



J d^y y 2 (P (r, 0) O(y , Vl ,y) P *(0, 0)) 
/ d^y(P (r,0)O(y , yi ,y)P^0,0)) 



(37) 



Based on (36), we obtain 



w 2 (r,f3, yi ) = -2(d-l)-^k,g 



M m M n 



y2M n + yi M„ 



V g - (3/2 M m + yi M n + 2^ ^2- 



2 M m y{+M n y 2 ^F m ' n (q 2 ) 



q=0 

(38) 



where yi = r — yi . At yi = 2/2 = 5 > the expression simplifies slightly, 
■■ 2 iY 1 ■-■ -■-){,!- 1) — 

where /i mn is the reduced mass of M m and M n defined by = M" 1 + M" 1 



-2(d - 1) 3^9 log { ^ b m b* n ^mn *V 



n ! F m ' n (q 2 )} , (39) 
J 0=0 



A. Effective string theory prediction 



On the other hand, Allais and Caselle [16] (see also the recent two- loop result [7J) obtained 
within the effective bosonic string theory the leading-order result, for x\ = 2y\ = r, 



d-1 r 1 Z 2 (f3,r) 

log — + — log » 

Z7TCT Tq 7TO" Zo(2p,r) 



(40) 
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Written in this form, it is clear that the second term can be interpreted as a difference of free 
energies. 



B. Transverse structure of the ground state of the closed string 



The limit y\ ^> (5 is most transparent from the closed string point of view, since the correlation 



function is then dominated by the closed-string ground state. Eq. (|38|) yields in that limit 
w 2 {r,P,yi) 



-2{d-l) 



yi (r- yi ) + (F°>°)'l 



2af3r ' F°>° 
where we have used the leading order relation M, 



-2{d-l) 



+ 



(F 



0,0y-i 



i?o,o 



(41) 



oL. The form factors are now evaluated 



at q 2 = 0, and the prime denotes differentiation with respect to q 2 . In the regime j/i > ft the 



effective string expression (40) behaves as 



lo g _^ + l_i r . + 0( e - 2 W/3 ) 



2tt<j 4r 4/3cr 



(42) 



It is consistent with the general expression (41) derived from the spectral representation of the 



correlator. The linear term turns out to agree automatically between the two expressions. From 
the closed-string point of view, this term is essentially a kinematic effect; we will return to its 
significance in the open-string interpretation of the three-point function. 

The rms radius of the closed string, defined in the standard way from the derivative of the form 
factor at the origin, can be identified with the r-independent term, 



("* ) closed 



2{d- 1) dF°>°(q 2 



q=0 



log . 

27r<7 4ro 



(43) 



F°.°(0) dq 2 

This term thus measures the logarithmic broadening of the closed string with its length j3. The 
prefactor is the same as for the open string, but the UV length scale appearing inside the logarithm 
is four times larger than in the open-string case. 

When r ^> (3, the open-string ensemble is at finite temperature 1//3. The local operator then 
probes the profile of the open-string states, averaged over with the Boltzmann weight. Equation 



(41) shows that the profile at mid-string grows linearly with the length of the open string [16\ 
This linear rise is likely due to the fact that 0(r//3) open-string states contribute to the correlation 
function when r and the width results from a stochastic superposition of these contributions. 

A linear increase is in fact nothing exotic, since for a screened potential V(r) ~ e~ mr , the profile 
goes like e -" 1 V( r / 2 ) 2+ 3 /2 , and hence the mean square radius is given by (d — for large r. 
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C. Interpretation of excited closed-string contributions 



Both the general expression (39) and the bosonic string formula (40) can be expanded in a series 



of exponentials that fall off increasingly fast. We require that the coefficients of these exponentials 
match. 

In the following, we use the leading-order relation between the matrix elements b n and the 
multiplicity (integer) factors w n , 

2 



l.o. Wn 
WQ ' 



(44) 



we choose them to be real and use the fact that wq = 1. We recall the values wi = d — 1 and 
w 2 = 1 + (d - 1) + \{d + l)(d - 2) 0, and also define AM n = M n - M . We start by analyzing 



the leading correction to (42), which comes solely from Zq(2/3), 



1XO 



(45) 



Expanding (36), one finds that the O(r) term cancels out automatically. From the 0(r°) term, we 

"Rei^V)' 



obtain the consistency condition 



d_ 
dq~ 2 



q=0 



2na ' 



(46) 



which dictates the strength of the off-diagonal matrix element between the lightest two string states 
at small momentum transfer. 

We now turn to the term of order e~^ MlT , which is of precisely the same order as e ~2^ M2r for 
the leading order spectrum. This time, both Zq(2/3) and Zq(/3) contribute and we find 



2wi -w 2 + \w\ 



2 "a AMir 



(47) 



while from the general expression, we extract 



w 2 {P,r,l) D -2(d-l)e~ AAhr 



dq 2 



F 



1,1 



W\ 



po,o 



+ 2JW2~- 



- ReF 2 ' 



2w] 



1,0 



ReF 
7?o,o 



(48) 



J q=0 



The comparison of Eq. (47) and (48) yields predictions for the form factor at small momentum 
transfer. By generalizing w 2 to values of y\ ^ U2, one could disentangle F 2 ' from F ' and obtain 
separate predictions for these form factors. In this way, a sequence of predictions are obtained for 
the form factors between low-lying states. 
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D. Three-point function with a non-scalar probe operator 



We now come back to (32) in the case of an operator with a more complicated tensor structure. 



Consider the case of O = T%\. Recall that direction '1' plays the role of time from the point of 



view of the closed strings. Then we replace Eq. (30) by (see Eq. ([18])) 



(m,p'\0\n,p) = \[E m (p') + £7„(p)] 2 /rV) " V/TV) (49) 
In this case we can write 

(Po( Xl ,x)0(yo,yi,y)PZ(0,0)) = £ ^ J ^* e -*» ■ (50) 

[|/rV)(4 + ir 2 ? - 2q 2 fT' n (Q 2 )} Zmn(yi,y2,X,q). 
where y% is set to x\ — y\ at the end. We now note that at leading order for large r = 2y\ = 2^2, 
(Wi + 4- 2 ) 2 ^n(yi,y2, 0, q) ~ (M m + M n ) 2 X mn (§, §, 0, q). (51) 
Hence to leading order the preceding analysis still applies, with the substitution 

F m > n (q 2 ) /r n (Q 2 ) " J^jT'^q 2 )- (52) 
A special feature of the operator T\% is that 

A m ' n (0) = S mn , (53) 



since the states \n,p) are energy eigenstates. In particular, Eq. (46) simplifies slightly to 



-2(d-l)Re(/ 1 1 '°)'(0) = ^. (54) 



It is interesting that, with our normalization of states ( 14 ) , the transition form factor is independent 
of p. 

V. LATTICE DEFINITION OF THE ENERGY-MOMENTUM TENSOR IN 

(d + l)-DIMENSIONS 

In this section, we derive the lattice form of the energy-momentum tensor in (d+ l)-dimensional 
SU(A^) gauge theory. Our main motivation is that these operators have been mostly studied in 
the d = 3 case, but recently there has been extensive work on strings in d = 2 SU(iV) gauge 
theories [T7j . This preparatory work will help us interpret those results. 
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We will follow the treatment [18J and generalize it to d dimensions. The idea is to identify the 
operators whose expectation value yield the thermodynamic energy density and pressure. We start 
from the Wilson action [T] on an anisotropic lattice |19| . 

S g = Y J PcS a {x) + P T S T {x). (55) 

X 

The action has two bare parameters, p a and f3 T , and there are two 'renormalized' parameters, 
the spatial lattice spacing a a and the renormalized anisotropy £ = a a /a T . At the isotropic point, 
Pa = Pr = P ( n ot to be confused with the symbol used for the closed-string length in the previous 
sections). The function S a and S T of the link variables Uu(x) contain exclusively spatial and 
temporal Wilson loops, respectively. The partition function Z depends on f3 a , f3 T and the lattice 
dimensions, N T ■ N d . The latter are related to its physical size by L = N a a a , Lq = 1/T = N T a T . 
We define the renormalized quantity Z by 

— N r 

log Z{p ot p T , N a , N T ) = log Z(J3 a , N a , N T ) - log Z(p o , p T , N , Nf). (56) 

The conditions that Z does not depend on a a or on the anisotropy £ translate respectively into 

LdlogZ L dlogZ ^ dp a d(3 T 
oL oL a log a a o log a a 

LpdlogZ \ - ^ , Q , <9/3 T /c , ,- a s 

where it is understood that the expectation values of S a and S T on the A^ ef • lattice are 
subtracted. We then recall the thermodynamic definitions of energy density and pressure, 

1 LpdlogZ 1 LdlogZ 

6 L L d 8L ' P dL L d 8L ' { 1 

With these definitions, we obtain at the isotropic point £ = 1, 

a d+ \e - dp) = J^{S„) + ^-{3r) t (60) 
d log Oi{j a log a a 

d d+1( \ ( dp a 1 dp a \ f dp T 1 dp T \ /oWfi1 x 

On the other hand, from the definition of the stress-energy tensor, we expect that 

(9) = = e - dp, (T 00 ) = e. (62) 

We also define 



1 



' flV 



T^u ~ j—^ S^6, (63) 
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so that in particular 



<fa> = ^i(e + P). (64) 



Since Eq. (60) and (61) hold at every temperature, we infer that 



a d+1 = J^S a + J^S T , (65) 
d log a a d log a a 

a + 6*00 = - — r + ttt, S ff - — — - + — — S T . (66 

\alog£ a + 1 a log a CT / yalog^ d + ldloga a J 

Recall that the magnetic field has d ^ d ~ 1 ^ components, while the electric field has d components. 

The lattice action can be expressed in terms of these fields, 

S(T = ^-Tr{B 2 }, St = ^Tt{E 2 }. (67) 

An important observation is now that at the isotropic point £ = 1, the operators 

S T - j^jSa and S T + S a (68) 

belong to irreducible representations of the cubic group in (d + 1) dimensions (32]. Since in both 
cases there is no other gauge-invariant operator of dimension [d + 1) in the same representation, 
both of them renormalize multiplicatively. 



A. The case d = 2 

Since the d = 3 case is well known |18|I19]. we focus here on the d = 2 case. The d = 2 theory is 
super-renormalizable, which leads to considerable simplifications. At treelevel on the anisotropic 
lattice, we have the following expressions for the bare parameters in terms of the renormalized 
ones, 

Ax = ^4, /3r = ^e (treelevel). (69) 

Hence 

W ' ^r— A, (70) 



d log a a d log a a 

8A - A- ^fir- (TD 



<91og£ " dlog£ 



Inserting these expressions into Eq. ( 66 ) , we get the following treelevel expressions at the isotropic 
point, 

a 3 9 = -p(S a + S T ), (72) 
a 3 9 00 = ^/3(2S C - S T ) . (73) 
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Since we already know that these linear combinations renormalize multiplicatively (see the remarks 
at the end of the last section), the full expressions for 9 and #00 read 

d(3 



<S a + S T ) 



'00 



dlog a 

lf3Z(f3)(2S a -S T ) 



(74) 
(75) 



with Z of the form Z(/3) = 1 + 0(/3" 1 ) and 



-/3(l + 0(/3 Now comparing these 



expressions with Eq. (|65j) and (66), we obtain at £ = 1 the relations 



<91og£ 

d((3 T - f3 a ) 
<91og£ 



d/3 
d log a ' 



(76) 
(77) 



Combining (74) and (75), 



a 3 T 00 



4 1 d8 

3 3 a log a 



)+Sr{ r \pZ(P) + \j£P-) 

3 3 a log a 



(78) 



By Euclidean symmetry, one then obtains also the expression for the diagonal stress operator, 



(79) 



and similarly for T yy . 

In summary, we have derived the lattice expressions for the renormalized diagonal components 
of the energy-momentum tensor. A simplification of the d = 2 case over the usual d = 3 case is 
that the one-loop quantum corrections to Z and d ^ amount to 0(a) effects, and the two-loop 
effects would amount to 0(a 2 ) corrections. The latter are parametrically of the same order as the 
usual 0(a 2 ) cutoff effects that are expected to occur in lattice gauge theory. For that reason, a 
one-loop computation is sufficient to yield a fully renormalized energy-momentum tensor. 



1. Application: width of the confining string 

In [5], the width of the string, stretched between two static charges separated by a distance 
r along the x direction, was extracted from the measurement of the Pq x = —So x + est plaquette 
expectation value at the midpoint x = r/2 (we now specialize to the case of the Wilson action; 
the additive constant drops out when subtracting the vev of the plaquette). We now interpret this 
result in terms of the energy- momentum tensor derived above. 
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Working at treelevel, 

(t \ 
1 -too > 

T 

XX 

\ Tyy J 

Inverting the matrix, one finds that 



+ 1 



-1 -l\ / S-A 
-1 +1 



V-i +i 



■1 



(80) 



5 



Sox 

Soy 



2/3 



2/3 



(Too + T x 



Now we use the general sum rules 



2/3 

\ j d d xT 00 (x)\^) 
/ d d xT xx (x)\iP) 



(Too + T t 



yyj- 



E. 



Lr 



dE 



(81) 
(82) 
(83) 

(84) 
(85) 



(V#> dL x 

Here L x represents an external parameter that E depends on. Thus, for a string of length r along 
the x-direction, 



(VIV) 



2 Q r 



<9r ' 



aE. 



(86) 
(87) 
(88) 



These can be viewed as the d = 2 version of the Michael sum rules |20j . For a long string, where 
E oc r, we expect the various plaquettes (summed over a time-slice) to come in the fractions 



(Sxy) '■ (Sox) '■ (Soy) 



1 



2' 



(89) 



We finish with a numerical application of Eq. (87) based on the data of For the rest of this 
section, we set the lattice spacing to unity. The profile obtained in [3] from the Sq x operator is to 
a good approximation Gaussian, with 



j — < 



dyAexp(-ly 2 /R 2 



2tt-A-R. 



(90) 



From Fig. (2) of [5), one reads off R « V12.1, and from Fig. (1), A « 0.00038. Thus the left 



hand side of Eq. (87) roughly amounts to 



2vr/3ARr. 



(91) 
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If we neglect the quark self-energies and the string corrections, E ~ ar and the RHS amounts to 

ar. (92) 

Numerically [5J, after simplifying the common factor r, we have LHS ~ 0.030 and RHS ~ 0.026. 
Given the approximations we have made, in particular the neglect of the quark self-energies and 
the use of the treelevel renormalization factors for the plaquette, the agreement is satisfactory. 

Based on the remarks of section (pd|), we expect all three plaquettes to yield the same string 
profile, to leading order in the string fluctuations: each of them contains a piece of either the energy 
density Too or the longitudinal stress T xx . This is indeed what the authors of [3] observed. We 



would however expect T yy , whose expressions in terms of plaquettes can be read off Eq. (80), to 
yield a different profile. 

B. The case d = 3 



It was numerically observed a long time ago [21J in d = 3 dimensions and for the gauge group 
SU(2) that the trace anomaly operator T w yields a large string width that grows with r in the 



observable (37). The linear combination 3Too — T\\ — (T22 + T33) was found to yield a smaller 
value, and no clear evidence for a growing width could be seen. According to our conjecture for 
the coupling of the stress-energy tensor to the string, this operator should measure the same width 
at leading order for a long enough string, but in the range ryfa < 2 reached in the study the 
corrections could be significant. 

A few years later, a new numerical study was carried out j221!23j in the same theory. The authors 
considered the operators T w and (Too + ^11 — (^22 + ^33) )> which yield similar profiles, as expected 
from the leading terms in Eq. ^ and Q. The operator that the authors call the 'transverse energy' 
is proportional to Too — Tn, and there is some evidence that the profile measured with this linear 
combination is indeed different, as we would expect based on the arguments of section (|TT|). The 
reader is reminded that we are using Euclidean conventions here, see the comment at the end of 
section ([lT|. 

VI. CONCLUSION 



In this paper we have analyzed the transverse structure of the confining string in non-Abelian 
gauge theories. By postulating a particular coupling of the stress-energy tensor to the string 
degrees of freedom, we obtained a prescription for how to compare the string profile obtained 
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from local gauge theory operators with worldsheet observables. We then derived the closed string 
representation of the three-point function from which the string profile can be extracted. For this 
purpose we first enumerated the gravitational form factors that characterize the string profile. The 
functional form of the leading-order prediction for the string's square width is then found to be in 
agreement with the closed-string spectral representation. Most importantly, we showed that the 
square radius of the ground-state closed string, defined from the slope of its form factor, grows 
logarithmically with the length of the string, just as the square radius of the open string does. We 
also obtained a prediction for the transition form factor between the ground and the first excited 
state. 

A number of open questions remain. It is not quite clear yet what role the ambiguity in 
the form of the worldsheet energy- momentum tensor (canonical vs. improved) plays beyond the 
quadratic order. It is not clear either to the author how exactly a generic local gauge-invariant 
operator couples to the worldsheet fields. Finally, we wish to comment on the prospects of fully 
characterizing the QCD string's structure. In the analysis of hadron structure, Generalized Parton 
Distributions have provided a powerful way to characterize the structure of a relativistic bound 
state such as the proton (see |24j for a review of the subject). Their moments in the longitudinal 
momentum fraction are given by the form factors of twist-two operators and are thus computable 
in the Euclidean theory |25j . These moments can be interpreted as Fourier transforms of the 
transverse distribution of partons |14j . The higher the dimension of the operator, the higher the 
longitudinal momentum fraction of the partons that it is measuring the transverse distribution of. 

It is a fascinating program to think about an analogous comprehensive way of characterizing 
the structure of the confining string. Here there is no need to go to the infinite-momentum frame, 
which leads to kinematic simplifications in the proton case, because the string is parametrically 
heavy compared to its transverse width. This warrants the interpretation of form factors as the 
Fourier transforms of 'parton' densities. By analogy with the analysis of proton structure |26j . 
a rational for which tower of operators to concentrate on may be provided by a 'Deeply Virtual 
Graviton Scattering' gedankenexperiment. Higher-dimensional operators will presumably corre- 
spond to probing the transverse distribution of 'gravipartons' which carry a higher fraction of the 
string's energy. By the well-known arguments, we would expect to find smaller transverse radii for 
these operators. 
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